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The  restfaroh  reported  in  this  doomont  has  been  mode  possible  through  the 
support  end  sponsorship  of  the  U.S*  Department  of  the  Army.  throv\gh  its 
Huropeon  Research  Offioe.  This  Report,  not  necoswailly  in  final  fioxm,  is 
lntcnde<l  for  the  internal  monagoraent  use  of  the  Contractor  and  the  U,S» 
Deportment  of  the  .Army. 


1.  QOffiRAL  SURTEI 


During  the  period  lot  Ifoveoibor,  1961  to  31at  Ootoher,  1962,  work 
has  continued  on  the  follcnring  projects: 
l.l)  General  Manv-Eleotron  Theory 

The  method  of  generalized  product  f^motlone,  developed  in  previous 
Reports  and  p\d>llcatione,  has  been  ertendod  in  two  directions,.  First,  the 
strong  orthogonality  condition  on  the  funoticns  describing  different  electron 
groupa  has  been  relaxed,  in  the  two-groiq)  oaso,  and  the  oona^iquertoes  of 
non-orthogonality  are  being  investigoted  in  detail.  Secondly,  been 

shown  that  when  one  of  the  two  grovp  functiems  belongs  to  a  finite  aubspooe 
of  H3,lb«3rt  space  (e.g,  an  antisymnetrized  orbitaj.  product),  the  results  easily 
reduce  to  those  of  Sinanoglu  and  Saass.  Some  attention  is  now  being  given 
to  the  determination  of  optinun  pair  functions  to  describe  pair  ooirelatioiis., 
Oji  the  other  tmd,  the  general  oaae  of  two  norr-orl'hogonal  gi*oi?)3  is 
sufficiently  ln^ortant  to  merit  tVu'ther  development:  it  is  of  special 
iJuportciTice  T^hen  the  two  electron  groups  lA  and  B,  oay)  ovo  ir.  aonrsinglot 
states,  for  there  is  then  a  non-trivial  spiiv coupling  problem  and  this  lies 
at  the  root  of  theories  of  ferrotnagnetic  and  anti- ferromagnetic  coiqplingj 
Considerable  progr  ess  has  been  made  in  this  area  and  it  ha.s  been  a]x>v/n  that 
the  Heisenberg  coupling  is  the  first  term  of  a  spin  Hamiltonian  repi-esenting 
the  interaction  energy  in  powers  of  the  spin  scalar  prx>duct  S'*^.  S®. 

Actual  applications  of  generalized  produst  Ibnotions  are  discussed 


in  other  aeotions^ 
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3.2)  CaloulationB  on  Atoms 

M  + 

fforic  on  the  8eq.uenoe  of  berylliun-like  systoaa,  li  j  Be,  B 
b8L8  been  concluded  and  a  paper  has  been  svtoitted  for  piiblication.  ?he 
calculations  suggest  that  the  electronic  stnucture  of  atoms  can  be  well 
repreaonted  by  a  wave  function  of  the  form: 

where  the  various  factors  describe  the  K,  L,  M, shells  and  raay  individ¬ 
ually  admit  electron  oorTfelation.  This  v?ork  is  fully  described  in  Special 
Technical  Report  SIR  No.  1  (1st  Sept.  1.962). 

1.3)  Caloulationa  on  Moleoulea 

Vtoi'k  on  the  NH^  radical  has  been  concluded  und  papers  are  in 
preparation  for  pub3,ication:  tho  work  has  alixsady  been  desoiubed  in  deta;Ll 
in  STR  No,  2  (15th  Sept.  1962).  The  calculations  have  been  concerned 
mainly  vrith  the  SixLn  density,  which  determines  tho  hyperfine  coupling 
constants,  and  a  variety  of  methods  have  been  eriroloyed.  These  include  the 
unrestricted  Hartree-fbck  metliod  (with  and  without  spin  projection  to  a 
douiblefc  state)  and  the  configuration  interaction  oet'nod  referred  to  in 
the  Pinal  Toohnical  Report  (30th  Nov.  I96I)  on  the  previous  contract, 

Non-eapirical  calculations  on  soall  organic  noleculcs  (e.g.  fbrtsal- 
ddyde)  are  being  continued,  vTith  a  view  to  developiaent  of  satisfactory 
semi- empirical  methods  ibr  dealing  ^vith  niuch  larger  molecules,  both 
saturated  and  unsaturated.  Tliis  involves  the  detailed  treatment  of  the 
<f  electrons,  and  o  group- function  approach  (each  localized  d  bond 
being  a  "^electron  garoup")  seans  to  offer  considerable  promise.  Owing  to 
tT>c  air  ties  ir  the  vcluos  of  many-cento’e  2-electron  integrals,  it  has 
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boon  neocssoiy  to  start  the  systsovxtlo  devoloincnt  of  inte^al  prasraniaes 
foi-  ©IsotrMilo  ooaputation*  Progreea  has  consequent iy  been  but  when 

all  prograaiaoB  ora  operating  it  should  be  possible  tc  deal  veiy  ragpidly  with 
a  Timber  of  Intereatijig  moleoulco. 

An  investigation  of  otoric  effects  in  anall  nolcouloa  is  also 
dependent  on  the  availability  of  Integral  values.  The  oaloulationa  have 
however  been  pressed  to  a  point  at  whioh  the  twiin  faotoj’s  deteroinlng 
inters  and  intra-moleoular  repulsions  ax-e  becoming  evident.  In  this 
work  a  sii^plified  oodei,  consisting  of  two  hydi'ogen  molecules  (repi-esonting 
two  nelj^bbouriiig  bonds)  is  being  considared  in  great  detail, 

Finally,  work  on  aeproxiraate  self- consist ent  field  theory,  as 
fi^jplied  to  the  ir  oleotrons  of  oonji^tcd  cyateoD,  has  continued.  One 
notable  advance  has  been  the  development  of  perturbation  nethodo  which 
preserve  sclf-conaiotsncy  r®  to  first  order  in  the  wave  ilmotion  (third 
ordoi’  in  the  energy) .  7/ith  auoh  oethodo,  it  is  posoiblo  to  define  atom 
and  bond  palariaablUties,  analogous  to  those  used  in  Ilflckel  thociy,  v/hjob 
permit  the  discussion  of  substitution  effects  \7itlJ.n  the  iVanowork  of  a 
self-oonsiatont,  mat^y- electron  thoorj', 

1.4)  Klectrpnlc  Coraputation 

DuriiVi  the  period  under  review,  much  time  too  been  devoted  to 
prograane  developmont.  iProgroBaes  new  include  integral  oooputation  (so 
far  up  to  three  oentroo  only),  natrix  trensibimtions  describing  change  of 
basis  ate,,  self-oonolstent  field  calculations,  and  various  steepest  descent 
and  matrix  techniques.  The  "Moroury'*  computer  used  so  fai*  will  shortly 
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be  r^laood  by  :the  "Atlas"  and  the  whole  oomputat tonal  situation  has  had 
to  be  reviewed.  The  Quontua  Vhoory  ftroup  is  to  obtain  its  ewn  anall 
oon^uter  (IBU  1^20) «  os  a  result  of  generous  si^port  from  the  NATO  Rescaroh 
Grants  PrograoDiey  ond  arrangemonts  arc  being  made  fbr  use  of  a  lauoh  larger 
machine  (Eivaliah  L'loctrio  "KDF  9") ,  close  to  this  University,  vAien 
nocessaiy.  Consequently,  progrooiaing  vtIU  prcjscntly  be  svrltohcd  Urom 
"Autexjode"  to  "itortran"  and  cxsnsiderable  redevelopment  and  rccrganisatlon 
v/ill  be  naoeasaiy. 

Vfork  Concluded  or  r-ublishoJ  (lat  Nov..  lS)6l  -  51st  Oct.,  1962) 

1)  Ih.nc  Currents  in  Arotaatio  Molecules.  I.  J.  Haddo:;  and  E.  MeV/eeny, 

J.  Cheu.  .  ^6.  2355, 

2)  Pertmtoation  Tlioory  ibr  the  S3Cl>Dirao  Density  Uatrix,  R.  ’.cV/eeny, 

Phys.  Rev. .  226,  1020. 

3)  Quantura  Organic  Cheaistry,  I?.  Ifc'.'eeny,  /uin.  Rep,  Chen,  Ooc.^  (icnrl.)  196.1, 
Z;)  Bond- length  Altemation  in  Heterocyclic  Rj.ng  Oystens,  D.  ,7.  Davieo, 

Nature,  194,  02. 

p)  Conforniaaione  Molcculare  cli  Poliraei’i  Idneari  in  Soluruone,  V.  liagnasco. 

La  Scuala  in  Aaiono  (Pcrioriical  of  Ihte  Najdonalc  Idi'ooarburi,  Italy) . 

6)  An  Empirical  Method  ibr  Calculating  Baiviox-s  to  Internal  Itotation  in 
Staple  i-Joleculcs.  V,  Liagnasco,  Nuovo  Cinonto  425* 

7)  7^Poly( ethylaie  oxide)  in  Solution,  (C.  itossi)  and  V.  Magnosco.,  J.  Polvacr 

977,  (1962)  . 

8)  Tha  Donsity  Matrix  in  Mary- Electron  Quantum  Mcdianios  Id.  III. 
Applications  to  the  bcrylliun  atom  euid  4- electron  iona.  R,  :bi7ceny 
and  T..  Sutolifi'e  (Suboitted  ibr  (-roc.  Hop.  Soc,,). 
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9)  Spin  Polarization  Effeots  in  Poraaognetic  i&^lcculos,  A  Colaulaticxi  on 
the  M2  Badical.  R.  McV/eoqy  £md  B.  T,  Sutcliffe,  Mol,  Pl^sioe  (to 
be  publiahed)  (Se©  also  STR  Ho.  2,  15t;.  Sept,,  1962). 
do)  Uhreztrioted  Hartre^-Jbck  Caloulations  on  the  KH^  Radical,  B.  T.  Sutclif: 
(to  be  published)  (Sco  also  STH  Hb.  2,  19th  Sept.  1962). 

11)  /stai'k  Effect  in  I^rogen,  E,  Steiner,  J,  Cliaa.  Fhys.,  45i. 

12)  T^Claaoical  and  Continuun  Contributions  to  rolorieobllity  and  Rtroe  ' 
Constanta,  (r/,  l^ra  Bi-own  and)  E.  Steiner,  J.  Chai.  Phye., 

2.  SCiamfTC  V/ORK 

In  this  aootion  we  give  flill  abstracts  01’  x'eccntly  published  papers 
and  Special  Technical  Reports  relating  to  woric  concluded  during  the  peidLod 
Ist  November,  1962  -  3lBt  October,  1962,  olong  v/ith  more  detailed  accounts 
cf  ^rk  still  in  ijroeresa.  The  v^ork  in  progress  rdll  'a©  continued  during 
the  cont.vact  period  oameocing  1st  Novoubor,  1962.-, 

2.1)  Wave  Functions  for  Borylllura  and  j-Electron  Ions 

The  group- function  lacthod,  and  the  forras  of  tho  inner- she! i  and 
outez^shell  functions,  wero  suraraarized  in  tho  final  Tcchnictil  Rcpoii 
(,30th  Novo.iber,  I96I)  on  the  previous  oontrcot^  The  whole  investigation  Isas 
nov:  been  described  in  detail  in  STii  Kb.  1  (Ist  Sept -jiiber,  I962)  and  a 
paper  has  been  subaitted  for  publication.  At  this  point,  therefore;  v/e 
cierely  reoall  the  results  for  the  boiylliua  atoa,  these  being  tjorical. 


Not  peribnaed  on  contract 
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The  oaXaulated  total  onorgiea  (in  liortre©  atcxnio  wits),  \7ith  and  without 
the  reotriotion  of  strong  orthogonality  between  inner  and  outer  ahells,  are 
oota£>ared  belov;  with  the  results  of  the  one-deteminont  calculation  (Slater 
orbitals)  vrith  Hartreo-Rxalc  culculationa,  and  v/ith  tlic  virtually  exact 
values  obtained  frm  exceedingly  large  variational  calculations. 


1  Derermlnant 

ihl' 

Or-tlcgonol 

Groins 

Non-Ort  logonal 
Groups 

» .  .  - 

licact 

-14.-556 

-14.575 

611 

-1k,621 

'■14*667 

It  is  soai  that  jnsooval  of  the  strong  orthogonality  requii'craerrt  leads  to  no 
groat  iraproverjcnt  of  the  grotq^  function  r3su3.t8,  cJ.tliough  it  introduces 
trcjaendoua  uoiaplicationa  into  the  calculation-  TIds  is  encoui’aging  raid 
ler-ds  sup-'-ort  to  the  belief  that  systemo  vd.th  scpaiv.te  pirts  can  ?)e 
aoou.ratoly  described  by  using  on  antis^-innctriced  prcdu-ct  of  group  iVnctionc, 

The  variational  calculation,  in  which  each  groviu  is  morle  sel-f-consiatent  in 
the  lisld  produced  by  all  other  groui:>s,  was  found  to  be  Ixjth  convenient 
and  ra^adly  convci^ent.  The  olccti'on  coxa-elation  introduced  within  each 
group  v/as  investigated  by  calculating  the  :xxir  function  for  elcotrcne  of 
opposite  Sidn  and  cettparing  tlds  vdth  the  (unco’/rclatod)  product  UiipixDxiitiation. 

The  resultant  "correlation  Iwlcs"  (radial  in  the  inner  shell,  angular  in 
the  outer  shell)  v/erc  plotted  and  discussed. 
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=  OOB©  + 

=  sinfl  - 


sin  9 

008  0 


(bonding)  (i  ^  1^2) 


(antibaadin#^ 


(4.2) 


The  ‘ground  state’  of  each  bond  is  then  described  by 

Bi(l)B^(2)  @Jl,2)  (iw3) 

while  the  lowest  triplet  states  are 

[B^(1)B.(2)-B^(1)B^(2)3  0,Jl,2)  (b  =  -1,0,1)  (4.4) 

Here  ®4,„  is  the  spin  ftsaction  oorresponding  to  spin  eigenvalues  s,m. 

The  triplet  funotions  of  each  bond  may  be  ooupled  vd.th  the  odd  oleotron  to 

give  an  eaooited  dciblet  figaotion  fbr  the  whole  Dolcoule  and  this  mixes  with 

the  ono-cc«ifiguration  ground  state,  putting  q)in  density  into  the  bond.  This 

2), 

is  the  neolioniso  invoked  by  several  aufchora  (o.g,  McConnell  and  Chestnut 

5) 

whioh  has  been  treated  oleewhcie  in  a  vory  general  way., 

The  (Kperimen tally  observed  quantity  is  the  cot,?>ling  constant  J 

•!V 

wliioh  appears  in  the  part  of  the  spin  Ifaailtonian  linear’  in  the  nuclear 
spins,  namely 


Tlii.y  is  related  to  the  spin  density  at  nucleue  n  by 


jOKa 


where  (Jj^(l^)  denotes  the  diagonal,  elarcnt  (l*  *  l)  of  Qj^(ljlO;  evaluated 
at  ntxjleus  n.  If  the  elootronio  g  ftretor  is  token  to  be  2,0023  and  the 
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spin  density  is  meosurod  in  atonlo  tsilts  of  volume  this  beoomes 


Jn  -  400.12  g„  Ql(g 
4) 


Ifc/seo 


The  nuclear  g  vuluos  assumed  are 

(hydrogen)  «  5.585W>  (nitrogen)  =  0;i»0357 
and  Tdicn  these  mlues  are  used  tiie  results  listed  in  Table  1  arc  obtained. 


Table.  1  Couplinpt  Constants  (Cl) 


Slater  AO's 


SCP  AD '  B  Expcriiaenta]. 


Nitiog«i  22.3  ll-,5  14.45  ife/sen 

Hydrogen  -38.8  -22.8  (»)  33.52  iic/sco 

( ii)  Unroetriotsd  H’>jtree-tbck  Calcul-ations 

in  this  apprO-ximation  the  electrons  nre  divided  into  on  -aoixi 
group  and  n  /S-apln  gTOT.?i,  and  we  take 

$!(«.. -N)  =  ...  w)J 

where  ^  are  usucdly  antisyoroetrized  orbital  products.  Ilie  spin 

denoity  is  simply 

Qj^(l)  -  (1)  -  (1) 

whore,  for  eacaaplo,  (l)  denotos  the  diagonal  elcncnt  of  the  apinloas 
+  + 

density  matrix  Pj^Cl;!)  for  ©lectz'ons  of  sp'in  +^.  The  variational  tact  hods 
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fbr  detemlning  and  (without  orthogonality  reatJlctionB  on  the 
orbital  factors)  arc  well  known:  they  comprise  tlK;  "unrestricted  Hartrcc- 
Rx>k"  approach.  The  ooupling  constants  detcminod  in  this  way  are  given 
in  Tdblo  2. 


Table  2.  Oouplinit  Constants  ftffiP) 

Slater  ilO’ 8  SCP/iD^s  Experlaental 

Nitrogim  75«2  39.4  l<!t-.45  Mo/aeo 

Hydrogen  ~35.3  -25,0  (-)  35.52  Jfo/soc 

One  tiia;}or  dicadvontoge  of  this  approach  is  that  the  functions  used 
is  not  a  pure  doublet ,  and  it  is  not  yet  door  how  important  the  contamination 
by  fVinctione  of  lilghcr  multiplioity  may  be,  ^inothcr  arfcward  featxiro  is  the 
"splitting"  of  the  Is  Core  orbitals,  vAiioh  recognises  oloctron  correlation 
in  the  core  (decreasing  the  variational  energy)  at  the  sane  time  as  any 
"spin  wlarization". effects.  Strictly  the  tv/o  effects  should  bo  independent, 
and  some  of  the  core  splitting  (whioh  appears  cuj  the  difference  of  two  very 
large  quantities  ond  is  thei'ofore  difficult  to  calculate  accurately)  my  in 
foot  by  sirorious,  IVo  modifications  were  tl\ercforc  investigated:  first, 
the  "freozing"  of  tlic  core  («x  and  ^  electrons  being  given  the  same  Is 
orbital;  and,  soocndly,  projection  of  the  onoconfigxiration  function  to  a 
pure  doublet  state.  The  first  modification  rras  of  limited  iiaportanoe.  But 
the  Second  yielded  rather  different  results: 


-  ll- 


j(nitrogcn)  =  26.2  j(lViiix)gon)  =  -16.7  (SCF  AO's) 

The  spin  ocntaainntion  hy  higlier  multi pleto  is  tliorofoi'e  a  very  oexlous 
factor,  i'iiia  contamination  diould  of  course  roMvod  before  the 
variational  calculation  but  this  leads  to  many  difficulties.  One  gi-eat 
advantage  of  the  oonfiguration  interaction  method  is  that  these  diffioulties 
are  avoided  fixm  the  start* 


1)  R.  ifoffeeny  and  K,  Ohrio,  Proc,  Roy.  Soc.,  A  25^,  367}  (i960), 

2)  If.  M.  MoConnell  end  D.  B.  Chestnut.  J.  Chan.  Pliyo.,  107,  (1958)* 

3)  R.  MeV/eeny,  Preprint  No.  57,  Quantum  Chetnistry  Cre-up,  Uppsiila  University 
(lith  Nov.  i960). 

4)  J.  Poplc-.  \V,  Schneider  ond  ii.  Bernstein  "High  Resolution  Nuolatr 
Kagnetic  Jlecotiance".  (jfcGi-av.'  Hill.  1959) 
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2.3)  Grpup-^\motion  Culoulatlcios  on  acall  ifolooulea 
MoV/ocny  and  Olino^  have  shown  tliat  in  the  caae  ofthev/ater  nolconle  the 
results  of  L0M>1S0  SCF  oolaiJntlcns  cxin  easily  be  surpassed  in  ooouracy 
by  using  on  antisymotrieed  prsduot  of  group  fimotions  and  including  Cl 
vrithin  caoh  of  these  grotf>  fliziotions.  The  aim  of  tho  voiic  desoribed  here  is 
to  cBctaid  these  oalculations  to  other  aoall  moleoulcs  and  to  investigate  the 
possibility  of  a  seniecsplrical  developoent  which  would  ooke  the  method 
applioablo  to  larger  moleoules.. 

The  starting  point  for  these  c  edculations  is  a  basis  of  orthogonal- 
ized  atomic  orbitals,  so  ohoeon  os  to  emphasize  the  locLiliaation  of  '  bonds 
and  lone-pair  electrons.  Eie  total  wave  funotion  is  then  given  by 

(1, 2, . . .  .N)  =  Mi^(-l)^\Kj^(lr 2)  \|/g(3,4)  , . .  (1) 

where  the  tv/o- electron  iUnction  describing  a  looclizetl  o'  bond  (R)  is 
=  Ci§i(l.2)  +  C?§2(1,2)  +  C5^f(l,2) 

with 

’.There  =  2  =- 

«2(1)  =  z'^CXad)  - 

Here  stands  for  2-Sdet  i;^j^(l) 06(l)  ^j(2)/3(2)}  and  and 

are  orthogonal  .<^’3  on  atom  a  and  b  respectively. 


/  (2) 

j 
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?he  total  electxmlo  energy  can  then  be  written  as 

vbero  i«  the  total  electronic  energy  of  tho  aystcra  with  group  H 

rcnoved,  while  Bg  ie  tho  enerffi^  of  group  R  In  the  field  provided  by  the 
rest  of  the  anlooulc.  Eg  is  obtained  by  solving  a  3  x  5  secular  equation, 
tho  matrix  elonents  of  wliloh  oon  be  derived  in  the  usual  my  in  tcnaa  of 

and  (lllm)  =  /i^(l)^(2)  ~  , 

trtier*  h  ff(l)  =  h^(3)  +  21  (  ^  ^(l)  ). 

h^(l)  is  the  Idnetio  +  potential  onet‘gy  oixjrator,  t7  (l)  and  (l) 

tho  CouloQo  and  exchange  operator  respectively  (of.  MciTcony^^) ^  Now  I.  . 

and  (ijjkl)can  be  oxpreesod  in  tenns  of  the  oori'esponcUnc  quantities  (X<v»Pob 

and  (obied)  over  AO' a  etoa,  where  the  ex’s  and  nev/  contain 

N  N 

a3.1  interactions  between  different  groupc  and  end  t)ic  bare-- 

N  r  fi  N  K 

nucici  quantities  =  JXo^  X  a'*'"  "  7x1’^  X 

=  0!'^  +  ^  {[(aaiuu)  +  (aalw)  -^(uulau)  — ^(avluv)] 

S(¥^K) 

■i-Ya  [(ooJuu)  -  (oalw)  ~i(aulau)  +  -A-(avjav) 
C2(aalu-v)  (au(av)]j 

Pab  '■  ^  ([(abiuu)  +  rablw)  -ii(au)bu)  “-^(avlbv)] 

•J-Yg  [(abluu^  “  (ablw)  -^(aulbu)  +-^  (avlbv)  j 

+X^,  (2(Qb|uv)  -^(aulbv)  ~^(ay{bi^3j 
R  R 

a  and  b  denote  AO's  and  involved  in  bond  R,  u  end  v  the 

AO’s  and  Xv  involved  in  bond  S. 
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Uie  paimaetera  X  and  T  ore  giv«n  in  terns  of  the  coefficients  in  aqu.  (2) 

s  s 

ty 


Xg=  (0^  -  (0^^ 

Tg  -  ^ 

Yg  reparesents  the  nisount  of  diazss  tmnsfexrcd  fron  orbital  v  to  orbital  u 
bond  3,  vhcreaa  Xg  is  associated  r/ith  the  o^wuige  density  within  the 
bond.  This  is  evident  from  tho  one-body  density  matrix 

pJdjlO  =  2/$g(l,2)  $g(l‘,2)dx,, 
which  becomes,  after  integrating  over  the  spin  variables  and  introducing 
the  quantities  d^ 


=■•  (1  ♦  Y^)  a^(i;iO  *  (1  -  dd') 

*  Xg 

Y,,  is  escpcoted  to  bo  small  for  bonds  of  lo?/  polarity,  being  datciinined 

essentially  by  the  difference  of  the  electronegativities  of  the  atoms 

to  wbich  u  and  v  arc  attached,  whereas  1X^1  is  expect td  to  bo  close  to 

b 

unity  for  looolisod  bowis.  (In  fact,  preliminary  calculations  for  Ii2G=0 
shew  that  lx  1=  0,97  ond  Y^  =  0,08  fbr  the  CK  bond{  con’csixjntling  values 
for  the  00  ^  bond  being  0.94  and  0.19,  and  Ibr  t!io  or  bond  0.90  end  0.08) , 

This  analysis  can  easily  bo  extended  to  include  conjugated  molecules, 

i,e,  raoleoules  rMch  contain  one  grot^i  of  more  than  t\:o  delocalized  electrons 

(it  eleotrons) .  If  wo  describe  the  or  electron  group  Q  aa  usual,  in  tr  election 

theories  by  a  ono-deteminont  wave  I'unotion  the  ta's  and  in  the  cnci^ 

expression  for  this  g  xoLcp  arc  atill  given  by  etiuation  (3) .  This  shov/s 

that  the  u'a  and  jO’a  defined  in  (3)  are  cloaelj'  related  to  the  corrooponulag 
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qxicntities  in  the  oaai- empirical  theoi*ies  of  Pariser  and  and 

where  the  neglect  of  differential  overlap  and  the  opprosiraatione  involved  in 
evaluating  the  Coulotab  integraia  (aa|bb)  can  bo  Justiilod  by  a3Buniing  that 
the  boBlB  orbitals  arc  ortl>ogonalized  AO*  a  (of,  IMTCony^^)  ,  The  presoit 
formulation  therofbro  suggests  a  ncyr  cf  calculating  the  oc’s  ond  fi'a 
in  these  theories.  This  would  bo  of  interest  for  in  cstimting 

the  difference  betwear.  the  paronetera  for  the  nitrogen  ataa  in 


u  •  O 


Qnd 


In  the  yrosenoe  of  a  ir  eleotroo  group  Q  in  the  noloculo,  equation  (3) 
for  the  ot's  and  ^  '  b  ot  tho  tv70~eieotrc«  Qvovp  R  haa  to  bo  diangea 
into  0")  • 

H  {((aa(uu)  +  (oalw)  ~  'J(aulau)  -  §(avlav)3 

+Y„  [(aal  uu)  -  (oalw)  -  ^(auiau)  +  -g-(avlav)] 


a 


K 


S(-R.a) 


+Xp  [2(na\uv)  -  (aulav)]} 

kl 

'  [2(aalr3)  •*  (arias)] 


>^5 


ra 


(30 


uu 


/2>S, 


+  X-  {[(ablutO  +  (ablw)  -  -^(ouibu)  -  -^(avlbv)] 

S  (*  R,a) 

+Tg  [(abiuiO  -  (ablw)  -  .*-(aulbu)  +  g-(ovlbv)] 

+X  {2(abluv)  -  ^(aulbv)  -  -^(avibu)]} 

+  Zr^.  [2(ablra)  -  (arlba)l 
rj  rs 

Tho  R -matrix  is  given  in  tema  of  the  mtrix  T  of  the  u  orbital  ICAO~ 
ooefflcients  by  JR  ‘•TT^.  This  shows  how  -  by  treating  all  S’  bonds  and  lone 
pairs  as  loceliacd  two* electron  groUi*  -  tho  serai- empirical  thooriea  of 
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Poplc  &  t?ari3er  and  parr  oaa  be  extended  to  include  tlie  S'  electron 
frameirorh  explicitly  and  to  treat  saturated  ooinix>und£  as  veil.  One 
possible  \vey  would  bo  to  oaloulate  all  cne-clootron  integrals  rigorously 
(this  is  q}x±te  feasible,  because  the  cno-  and  tv/o-centre  integrals  oan 
be  cnroltmtcd  in  closed  form,  and  the  nunber  of  three*centr3  integrals 
is  rather  limited  and  to  get  the  oc  's  and  ^  's  (over  orthogonal  hybrid 
orbitals)  by  a  sytaaetricnl  Lflwdin  orthogonolization.  The  two- electron 
integrals,  which  are  very  nuaerous  and  the  evaluation  of  '.7hioh  is  very 
time  oonauning,  could  then  be  approxijaatod  alorig  the  lines  suggested  by 
l^i'iser  A  Parr  /Jind  Poplo,. 

The  Ibmaldehyde  molooule  \7a5  chosen  as  a  first  exanple  to  ooDipar'e 
various  appreatiaations  with  the  exact  group- fbnotion  calculations  outlined 
above,  largely  bocaivsc  the  integrals  over  Slater  orbitals  ax'c  available 
(Dunean  et  al.^^),  Progrotanea  have  been  written  6Uid  t  ested 

a)  to  generate  the  transformation  matrix  which  transforcis  the 
original  basis  into  orthogonolized  hybrid  oibitala  loca].isG<i 
in  the  bonds; 

b)  to  transfoini  the  ono- f^lection  integrals; 

c)  to  transform  the  tv/o~ electron  intqgreje. 

d)  to  perform  the  group- function  calculation  in  on  iterative  ?/ay 
by  solving  the  secular  equation  for-  aich  <  roup  in  turn,  starting 
■with  Y=0  and  X  =  1,  until  self-oonaistonuy  ie  roaohed. 

Several  modifications  of  progrorjae  d)  have  been  written  to  calculate  caccitod 

states  and  to  ti-eat  the  ir  eleotrcns  and  the  oxygen  lone-pair  eiijectrons 

together  as  a  four- electron  giX)Up  etc.  Pr-ogromnics  b)  and  c)  ai*e  completely 

general,  a)  and  d)  can  be  easily  modified  and  adapted  for  the  calculation  of 

other  aoleoirlcs,  Ko  final  results  Iiave  been  obtained  so  fhr,  errors  being 
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suBpeoted  in  aoae  of  the  integials  listed  by  Duncan  ct  ol.  (Sens  exoJ»nce 

integrals  beooac  negative  after  the  tranolbruation  to  on  ortioipnal  basis)  ^ 

Sane  ncra  tlirec- centre  integrals  are  aloe  required  in  tranafomiing  tlie  one- 

centre  integrals  to  an  orthogonal  basio,  Sheso  oancol  each  other  in  a  tjolooulor* 

orbital  treatment  and  are  therefore  not  given  by  Duncan.  In  jxocnt  months 

the  cnin  anphasie  has  therefore  been  placed  on  the  development  of  cou^^uter 

progronmos  for  tlic  evaluation  of  all  molcoulaa.  integrals  over  Slater 

orbitals.  So  for  progitannes  for  the  evaluation  of  the  follovTin.ij  intcgials 

have  been  written  and  are  in  working  oi'der; 

(i)  ono-centre  intcgmla 

(ii)  two-oentre  one- electron  integrals 

(iii)  two-ocntre  Coulomb  integrals 

(iv)  two- centre  liybrid  integrals 

(v)  tla^oo-cer.tre  one  electron  integrals 

Programme  (i)  xbllavs  the  anal^aia  by  .  prt)grai;i;.ies  'ii.)  and  ;,iii) 

that  of  Irogrcmct’  (iv)  and  (v)  use  tw.-)-  itnd  thrca-diracneional 

nufacrical  integration  (of,  Magnusson  and  .tauli^^),  TbsO  Ikimeti-Goulson 

10) 

method  is  intended  to  be  used  fo.r  tne  3’or.Tainir);.'  integrala,  as  v.'ell  ca 
for  the  threc-ccntrc  one- electron  integ.tals  where  it  is  toped  to  be  faster 
than  numerical  integration;  Gauss  integiution  using  C192  ijoint  lattice 
(16  points  for  both  variables  and  ,  32  for-  d,)j3ivcB  only  on  aocui'acy 
to  thi’ec  figures  for  a  throo-oentre  one  electron  Integra?,,  and  takes  approK- 
iraatcly  12  ain,  Iwo-diniensional  integration  using  a  net  of  512  points  give 
a  six  figure  aocurrxiy  in  the  ease  of  hybrid  integrals  and  takes  ap-proximatoly 


20  sec. 
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2.4)  Inter  and  Intra-!!olecular  Interaotlons 


a©  problen  of  tho  intoraotlon  of  two  Qg  aoleculoo  i;  being 
Inveetir'atod  as  v.  first  step  in  the  study  of  lore  oocoplionted  oyatens 
(e.g,  OTj  -  CHj  nnd  tho  related  problen  of  bcrriors  to  internal  rotations). 


Tlio  frco.’otr  of  the  ayste-.  is  illuotrated  in  Figuro  1  #  One  nolooule 

io  fizod,  the  o<.her  ono  can  rotate  around  the  ac  c:cis  laoninr:  the 

dlctanco  D  s  R  con  tnnt,  its  ?tp.te  of  internal  rotation  boinr 
ao 

snooifiod  an  ajigle  9  ncaaurod  from  tlie  coplanar  position  in  nliioh 
tiio  nneloi  b,d  l^e  on  tlie  oane  side  -..itli  rosnect  to  tlio  lino  oining 
the  nvcloi  a»c-  (  i  =  0*),  R  i”.  the  eqiiilibrinr.  internrclonr  dlotanoO 
of  a  rig  noleculo  in  its  rro’>nd  3ta.te  (^E^).  Our  nain  interest  in  in 
voluQB  of  R.  A"  a  basis  n©  use  4  atonic  orbitaln  a.h.c.f  which 
have  the  fom: 


J 


20 


i.e.  Ifl-bydrocen-lUsa  A0*s  oentred  on  the  nuoloi  a,b,c,d.  The  ayiabol 
1  atcnde  for  the  t-'urae  apatieJL  oo-ordlnntoo  of  electron  1,  la  tho 
dlstcAoe  of  eleotron  1  froo  ^  nuoleua  n,  r.nd  I:  Is  the  uauol  sereeninc 
oonatont.  ’.7e  rjro  atudyinc  trro  different  appror.o’MB  to  thin  probloat 
a)  The  onercy  ^.p^rocoh. 

B)  The  force  field  anproooh. 

a)  Oio  onerg/’  approp.oh 

’.7e  aosune  that  t!io  trro  ayatema  A  and  B  (each  oonaiatina'  of  tr;o  11  nuclol 
and  tr?o  electrons)  ai^  approxL'nr.toJ.y  oeparr.'ed,  in  the  oence  thnt  the 
apatial  parts  of  the  xmve  functions  dosorlbin,  tho  aeparc-tod  .aoloculea 
from  a  qur.ntvin  eohanical  point  of  vior  da  not  h  ve  an  apnaaeciable 

1)2) 

overlap.  Tlie  validity  of  tills  approxioatlon  Iir,3  been  disonoaod  elsenloro.  ' 

Lot  va  suppose  tliat  A  end  B  are  In  their  oleotonio  pround  state  ;  and 

5 

that  wo  can  describe  eooh  of  then  b  -  a  "Ileltler-London  +  ionic"  ’.vavo 
function  of  the  f'oznn; 

fg(3,4)=[Xifcc(i-c|c)+  \zCC 


tdiere 


are  indiylduolly  nornioliBed. 


’7e  conotruot  the 
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wave  function  ^  for  the  total  uj'oter.i  in  the  fonaj 

lf^(l,2,3,4)  ^  W  C^§  U,2-)fg(3,4) 


there  U  is  n  nornialieation  factor  and  ^  is  a  partial  nntisyTiotrizor, 
30  that  ^  is  coaplotely  cjitisynj.iut6r<c  t;ith  roopoct  to  the  oxclianco  of 
any  pair  of  olocti'OiTS.  tfo  oan  at  this  point  distincnio'i  bct/eon  trio 
ca-.eo: 

l)  rtronc  ‘^rt  ofonal  case. 

II)  Ilon-ortltofonal  case. 


l)  '’■tronn  ortioronal  case 

Vq  otart  froD  an  ortho-norrmllzed  oot  of  AO*b 


IP  -  [a  b  G  d  ) 


TThich  may  be  obtained  fron  tlic  non»orthoconal  sot  jp  =  (a  b  c  d)  by 
moans  of  a  Bynnetrical  ortlioconalization:^ 

■~^/z 


f-n 

nhere  %*PP  is  t!iG  4x4  overlap  matrix  in  the  basis  p  . 

Wo  conotruot  from  (a  b)  cmd  fro©  (c  5).  In  tliis  case  the 
ccnditiona  of  ntronc-ortiiofronality  holda; 


Th®  Ikuailtonion  for  the  systein  nay  bo  mritton  as: 

■4=(r  ~  *'■ 

a. 

i,=l  '< 

i(o=-iv"w^v,(i)-v,(o  =  ^r. 


v;hore  tho  ayTibols  Ir.vo  an  obviova  ixsanin;;  'nd  atomic  ntis  (an)  arc  nood 
throuchont,  Tl'jo  oiqiectation  value  for  tlie  olectronic  enorc/-  is  given 
A  B  A  B  A6  AB 

F-E+E  +  -K 

^  '“o  o 

whore:  .  *  cr  ^A 


^\Ci)  ii)i) 

jA®  Jf 

k"*! 4  jr^az)  P/(2;l)  P®ai2)  ^lao/ir^  . 


A  /  A 

liar®  ^  rospectlvoly  the  1-body  end 

the  2-body  sniiless  de  eity  natricea  for  the  ayoton  A, ^(1)1) , 


the  soB»  qiuintities  for  tlie  systen  B.  Every  tern  In  t!zo  preceding 
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egression  for  tlie  energy  hns  a  very  8ii!5>le  pliysical  neaning,  J'or 
compute tional  purposes  It  Is  convenient  to  write  the  preceding  equ-.tion 
in  a  -vitr'x  fom  as  follows: 

t-r  Bn  Pi  +  4  ^  ^  > 

whore  *tr*  'neano  "trace",  and  h  and  g  are  respectively  t!ae  mtrices 
collecting  the  1 -electron  and  2-eleotron  intogralo  in  tlio  basis  p, 

6) 

They  are  related  to  the  non-ortliogonal  basis  p  by  the  raatrix  rolations  ' : 

oil  flu 

Here  jig  is  the  4  x  4  natrix  of  the  i -electron  intogralo  in  the  basis 
p  ana  j^tho  10  X  10  2-elootron  nv-trlx  in  tho  sane  basin  (the  riatrix 
oloinonvs  bo'ng  iirranged  in  a  dictionary  order),  Hhile  this  aotliod  is 
r/ell  adapted  to  r.  calculation  of  the  total  electronic  energy  (end  easily 
extensible  t )  nore  oooplioated  systeas) ,  it  does  not  allor;  for  studying 
explioity  the  dopondenco  of  tlie  electronic  energy  on  the  rotation  angle  0, 

II)  Mon-ortliogoral  case 

In  this  case  v/e  use  an  a  basis  the  set: 

jp^(oo  b  C  d)  . 


Ihe  electronic  interaction  energy  is  given  by  j 
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2  ab  a 

■*■ -t^n  ^  Sh  (.fo’*'Eo^J 

8L  Aa  '''■"■■ 

E  5h 

1 


where  the  tonne  L^s  describe  tlio  r:o  called  "nonctration 

n  n'  tt 

effects"  duo  to  the  overlap  of  the  electron  clou's  ~  and 
represent  in  tliis  case  tho  onercies  ef  truly  ooparnted  s  stens  A  and  B, 
and  -ro  independent  on  0.  For  sufficiently  lore©  values  of  D,  tlie 
■'rocedinr  equation  :x?.y  bo  ’.vitli  a  ipod  a.pproxi- 'r.tion  I’eiresentod  b;r; 

/  A  ^6  ,.A&W,  r.Afe  AB  ^AB  .  8>n 

Ej  <w(C  *c  +7  -K 


Usinc  the  Hullilcon  appro::imation  for  evaluatinr  ■tI\o  molecular  intocrals, 
v.'O  dind  th't  the  total  interaction  onorcy  takes  the  fom: 

E  -  e'  +  4(r)  +  1/-^ 


3) 


where  r  a  ®bd  (au)  ~  j;’  and  f^  are  constant,  f^(r)  and  f^Cr)  are  fVnotions 
involving  pol’/nonials  in  r.  Because  tlie  variation  of  and  with 
r  is  nuoh  smaller  than  the  exponenxial  one,  we  can  asoi’no  f^(r)  and 
oonsthnt  with  respect  to  0  and  equal  to  an  average  value  f^ 

So  we  get  the  expression: 

E-e'+  -fi 


~h 


i 
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This  focprossion  Is  olosoGLy  ainllor  to  a  preoedlog  relatinuahlp  Ibund 

etaplrioallj'  etnd  vhloh  olloocd  to  c^doulate  the  haxrier  to  internal  rotation 

9) 

in  ethono-likc  Qoloouloa  with  a  vexy  high  degree  of  aocurao^/'  .  Calculations 
for  different  values  of  D  and  0  using  the  exact  eqmtion  are  in  progress 
for  testing  the  iciportanoe  of  different  approximate  wave  functions  and 
the  possibility  of  approximating  the  molooulor  integrals  by  the  llullikon 
appronioation®^  and  biinlar  eKiionsioa^^^  .  ibr  comparison  the  exact 
calculation  of  uultiocntre  integrals  by  the  available  standai-d  laethods^^^ 
is  in  progix'ss. 

B)  The  force  field  approach 

This  approach  is  based  on  a  theorom  due  to  liellmann^'^'^  and  Ibynman^'^,  wijioh 
states  that  the  force  on  a  given  nuoloua  in  a  laoleoular  syst-.a  is  Just  the 
olassioal  eleotroototic  Ibroo  tint  would  be  exerted  on  this  nucleus  by  the 
other  nuclei  and  by  the  electrons  regarded  as  a  static  oliaige  distribution 
of  density  Onoo  tlio  electron  charge  distribution  has  been  doter-iined 
according  to  the  mctliods  of  molecular  (luuntum  meohonios,  the  force  on  a 
given  nucleus  in  a  molecule  may  be  caloulated  on  a  purely  classical  electro¬ 
static  basis.  If  S  is  the  ttc^icctatlve  value  ibr  the  energy  of  a  syutom  with 
wave  function  ^  and  Hamiltonian  4^,  v/e  may  write; 


i 
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because  is  boraiteon.  Ibc  sooood  Intc^al  in  the  rl{;;ht-l)cind  side 
vanishes  Identically  for  a  stationary  state  uhun  ^  is  the  exact  ;mvo  iUnotion 
for  the  syatem  and  also  when  ^  is  a  variation  fVjootion  containing  suitable 

lc> 

poroaoters.  In  the  .ircsent  oase  v/e  uioyly  asaumo,  at  first,  ttot  tliia 
term  la  small  oomijared  v/ith  the  first. 

In  this  apjoioxiaa-tion,  v/c  get  tho  fjlla/ing  expression  for  the 
expectation  value  of  the  foroe  Fe  acting  on  the  nuol.us  d: 


booauoe  the  kinctio  t  n-rgy  part  of  II-  does  not  de.-end  on  tlK?  cngle  B 

An  inspootion  of  the  explicit  form  of  tlw  i-otential  cnci^  operator  V  shows 

-1 

tliat  the  only  toms  vlJloh  depend  U[>on  0  are  the  nuclear  repulsion  Ru 
bet^voon  tho  nuclei  b,  d,  and  the  nuclws- electron  attraction  ^ 


r 


.  ?7  •• 


! 


Lat  us  assume  a  syeteai  of  oo>oirdinates  as  that  in  FLcwe  2,  and  lot  us 
soparato  tho  foroe  In  the  part  depending  on  tho  nuolcu»*nuolous  intor> 
aotlon  and  the  part  d(^endlng  on  tho  nuolouo-elootroa  interaotion* 

\7e  get  oaall^i 


© 


The  last  caiu^jssion  can  be  easily  sioplifled  Introdxiolng  the  l»body  spinless 
density  matrix  P^(l;l)  for  the  uhole  systcxi: 

F 0^  sf  R  j  si>y^  9-  ‘ PiCl;l)citr;i 

The  eTa?.uation  of  the  expootatlon  value  for  is  reduced  to  tho  evaluation 

of  a  certain  nuobor  of  l-eleotroo  moleoulor  intograls  involving  1^;  and 

ihoentro,  and  this  ney  be  oonvenicntl,''  ^:)crfbmod  by  aeons  of  a  numerical 

12)  e 

integration  in  olliptioal  oo-ordinates,  '  Onoo  F^  has  been  evaluated  for 
given  values  of  ^  ,  the  energy  as  a  function  of  the  angle  0  oan  be 
obtained  by  integration* 
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'i7ork  alone  these  lines  is  in  ^sroeress  for  a  oaQx>arlscni  between  tho 
two  different  approoolicG. 


Beferowea 

1)  K.  Cr,  Parr,  P,  0.  Sllison  and  * .  G.  Lyiros,  JCu.’,  1106,  (1550  • 

2)  J,  U.  Parte  and  K.  G.  .^irr,  JCP.  23,  335,  (1953) . 

3)  1'.  0.  lA/din,  JCP,  1^  365,  (1950). 

4)  R.  MsWosny,  PES,  A  253.  242,  (1955) . 

5)  R.  MoWceny,  Rev.  :iod.  Pl^ys..  ^2,  335,  (i960). 

6)  R.  MoWeony  and  R.  Olino,  PRLi,  A  255.  36?,  (1560) . 

7)  R,  MoWeony,  Q.  Cheniciry  Groupf  Jan,  15th,  Prc^viint  No.  59,  Uppaala  1561. 

8)  K.  S,  :!Ullikon,  J.  Chin.  Phys.,  ^  500,  521,  (l'A9). 

5)  V,  j'ingnasoo,  Ihwvo  Cinento,  425,  (1962) , 

ID)  E.  A.  iiaaon  and  J.  0.  lixrsohfeldcr,  JCP^  173,  756,  (1957) . 

U)  Baiuett  M,  P.,  and  Coulscm  C.  A«,  Phil.  Truns.  Roy.  Soc.,  A  243.  221, 
(1951). 

12)  J,  A.  IlQgnusBon  and  C,  iiauli;.  Proo.  ?l^a,  Soc.,  53,  (1961). 

13)  E.  Hellrnann,  Quontonohoale,  Deutiohe  and  Co.,  licipzig  1937,  p.  ^5« 

14)  R.  P.  Psynmon,  Phys,  Revs,,  34O,  (1939). 

15)  A.  C.  Burley,  Proc.  Roy.  Soc.,  A  226.  17O,  179,  (1954). 


-  29  -> 


nmm. 

to 

FTR  on  Contract  I:o.  DA- 91- 591- -IXJ- 20 34#  0I-26235-B 
I>  Porawnol 

During  tlic  oontract  p^-riod  Dr^,  D.  Davica  hao  loft  the  Group  to 
take  an  appointment  as  Looturar  in  ihoorctical  Clier.iistry  at  the  Univoratty  of 
GrSnincen  (Nethcrlanda)  ,  Dr.  B.  T.  Sutcliffe  haa  joined  the  Group  at  II.I.1'. 
(uoricin"  v/ith  Frofeasor  ijamett)  j  and  Dr,  Aa^ar  Ali  hao  been  apix>intod 
i-rofoaaoi'  of  Physical  Cheniatry  at  ProBidoncy  College,  Calcutta.  Dr.  hiioh 
Steiner  hu.a  joined  the  Groijq?  and  it  is  hoped  that  further  appointment  a  \7ill 
b^.  made  U'-ulng  the  year. 

II.  uan  hours.  .Materials.  pro:)CX’l:v 

A  rough  oatiraatc  of  man  Ivoura  expended  on  the  project  during  the 
oontract  pciTlod  ia; 

2000  hrs 
1000  hra 
2000  hrs 

1000  rJTS 

2000  hra 

Prinolpal  Invest i^pitor  500  hra 

Only  one  of  tlicsc  apixsintucnta  una  provided  wholly  at  contract 
oKpenae,  others  being  maintained  wholly  or  in  ^art  from  other  acurccs 
(Fcllowaliips  etc.)  ^  No  expendable  raaterialo  were  involved  apart  £ix)U  those 
inoluded  in  ovorheeda  (stationary  etc.)  and  no  important  ^jroperty  v/as 
aoq,uirod  at  diroot  coot  root  cxpc3iao„ 


Dr.  A,  Ali 
Dr.  D.  W,  Davies 
Dr.  M,  Klossingcr 
Dr.  V.  hagnarico 
Dr.  A,  T,  Sutcliffe 


I 
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m  WOTHJ  PMIS 


Durisig  the  neact  oontxoot  period  let  November  1962  •  Jlst  Oo^ber, 

It  ia  h(^>ed  that  four  pieoee  of  noxfc  will  be  oonoludedt 
1}  The  invoatigatioB  of  the  intorootion  betveen  ooloouXoa,  in  a 

variety  of  aptjroxliiiatlonay  and  of  its  iaplicatlons  In  the  theoxy  of  barrier 

heists. 

2)  The  ootopletion  of  oomputer  progreancs  fbr  oaloulationa  on  amoll  moleouXea, 
vith  applioations  to  ibnaaldehyne. 

3)  The  derivation  of  tlie  apin  IlaDiltonlan  r;liioh  rcproae^ite  the  interootion  of 
tv/o  eleotronio  ayatemQ  in  the  poners  of  the  soulnr  product  of  their  spins, 

uf.  to  terms  in 

k)  The  oalculAtion  of  apinropriatc  self-oaneistoat  v>eleotrDn  polari::^abilitle0 
for  a  number  of  sliaple  oonjugatod  moloouloo. 

Long  tens  projeote  oonocxned  vritb  general  moiiyelootron  theory  and 
tho  oorrelatlon  problem,  and  also  edth  the  development  of  computer  progroanea, 
aill  of  oourao  oontinue. 


I 


